Dynamic effective mass of granular media 
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We report an experimental and theoretical investigation of the frequency-dependent effective 
mass, M(a;), of loose granular particles which occupy a rigid cavity to a given filling fraction, 
the remaining volume being air of differing humidities. The effective mass measurements allow 
us to study the mechanisms of elastic response and attenuation of acoustic modes in granular 
media. We demonstrate that this is a sensitive and direct way to measure those properties 
of the granular medium that are the cause of the changes in acoustic properties of structures 
containing grain-filled cavities. Specifically, we apply this understanding to the case of the 
flexural resonances of a rectangular bar with a grain- filled cavity within it. The dominant 
features of M(a;) are a sharp resonance and a broad background, which we analyze within the 
context of simple models. We find that: a) These systems may be understood in terms of 
a height- dependent and diameter-dependent effective sound speed 100 — 300 m/s) and an 
effective viscosity 5 x 10^ Poise), b) There is a dynamic Janssen effect in the sense that, at 
any frequency, and depending on the method of sample preparation, approximately one-half of 
the effective mass is borne by the side walls of the cavity and one-half by the bottom, c) By 
performing experiments under varying humidity conditions we conclude that, on a fundamental 
level, damping of acoustic modes is dominated by adsorbed films of water at grain-grain contacts 
in our experiments, not by global viscous dampening, d) There is a monotonically increasing effect 
of humidity on the dampening of the fundamental resonance within the granular medium which 
translates to a non-monotonic, but predictable, variation of dampening within the grain-loaded bar. 

PACS Numbers: 45.70.-n, 46.40.-f, 81.05.Rm 



I. INTRODUCTION 

Loose grains, made of a variety of different materials, 
dampen structure-borne acoustic signals very efficiently 
when they partially fill cavities within the structure itself 
[ll^SHl- ^^^^ ^his reason, there is a practical motivation to 
develop an effective method to optimize the dampening 
of unwanted structure-borne acoustic signals. The funda- 
mental origins of the dissipation mechanisms in granular 
materials are still unknown, however, making it difficult 
to optimize the effect. Partly, this is because, until now, 
it has not been possible to study the relevant proper- 
ties of the granular medium directly, independent of the 
structure whose acoustic properties are being modified 
by the granular medium. 

In this article we pursue the concept of the effective 
mass, M{(jo)^ of a loose granular aggregate contained 
within a rigid cavity [1, Q ; a preliminary report on some 
of our results has been published previously M{uj) is 
determined by simultaneously measuring the force, F{u;)^ 
on the cavity and its acceleration, a(cj), when the cavity 
undergoes oscillation at a frequency uo. If an identically 
filled cavity is located within an acoustically resonant 
structure, we demonstrate how to predict the changes in 
the acoustic properties of the structure, such as sound at- 
tenuation or resonance- frequency shift, based on a knowl- 
edge of M{u;). This fact allows us to focus our efforts on 
understanding the properties of the granular medium di- 
rectly, rather than having to infer its properties via its 
effects on the host structure. In this regard we demon- 
strate the dominant effect of humidity, both on M{u;) 



and on the acoustic resonance of a steel bar having a 
grain-filled cavity. 

There have been several previous investigations into 
the origins of particle damping. Cremer and Heckl Q 
concluded that damping is especially high when the (ver- 
tical) thickness of a granular layer is equal to an odd mul- 
tiple of a quarter wavelength in the granular medium. 
Sun et al [3 have treated the acoustic effect of the gran- 
ular medium as if it was due to radiative damping; they 
computed the acoustic loss due to radiation by assuming 
the granular medium is a low-velocity fluid. Bourinet 
and Le Houedec and Varanasi, et al. [1, 0] have 
each considered the acoustic propagation characteristics 
of long hollow tubes, partially filled with granular mate- 
rial. Each approximated the medium as a low velocity, 
high attenuation fluid and each achieved quite reason- 
able agreement between their computed and their mea- 
sured values. (The two theories are similar but differ in 
their details.) We do not dispute these aforementioned 
results which are, in fact, quite reasonable. We point 
out, though, that the relevant granular medium parame- 
ters (the sound speed, the loss factor) are generally set by 
requiring a match to the observed acoustic characteristics 
of the grain-loaded structure. (An exception is Ref.Q.) 
It is this last feature that we obviate in the present arti- 
cle: For the kinds of structures we consider here, acous- 
tic loss is determined by the imaginary component of 
the effective mass, M2{oo)^ evaluated at the propagation 
frequency (or resonant frequency, as the case may be). 
Moreover, we establish that the properties of granular 
media are very much dependent on the filling level in the 
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cavity and we demonstrate that the side wahs of the cav- 
ity hold up some of the dynamic load. Thus the granular 
material cannot be idealized, for acoustic purposes, as a 
fluid, and certainly not a homogeneous one. 

Intuitively one might expect that the effect of the gran- 
ular loading would be to lower the resonance frequency of 
the structure holding the grains. While this often hap- 
pens we show situations in which M(cj) is negative in 
the frequency range of interest, leading to an increase in 
the resonant frequency of the structure containing the 
grains. This behavior has been observed before by Kang 
et al [10] who monitored the resonances in a clamped 
plate as more and more grains were loaded on top of 
it. Initially the resonance frequency drops; it reaches 
a minimum, then increases, eventually often exceeding 
the original (unloaded) resonance frequency of the plate. 
This behavior has a simple understanding in terms of 
a resonance within the grains, whereby the real part of 
M{uj) can take on negative values. (See Section IVII B[ 
below.) 

Generally speaking, we flnd that M{uj) exhibits a sharp 
resonance, which, as one part of our analysis, we inter- 
pret in terms of an effective sound speed, albeit one which 
is dependent upon the depth of filling of the cavity, and 
a broad tail that decreases roughly as which we 

interpret in terms of an effective viscosity. These gen- 
eral features have been observed previously, by others 
d, [Bj] using similar measurements. It is to be emphasized 
that our experiments are all done in the regime of linear 
(small amplitude) acoustics; this viscosity is not relevant 
to a granular flow, e.g. in a pipe or in a loading hopper. 
Each of these two interpretations is based on toy models, 
which assume that the entire effective mass is borne by 
the bottom of the cavity or by the walls, respectively. A 
concrete example of the former behavior is provided by 
the effective mass of simple liquids; we demonstrate how 
our technique enables us to measure the density and the 
sound speed of four different liquids. 

We have developed molecular dynamic simulations to 
analyze the expected behavior of M{u;) under the as- 
sumption that the contacts are described by dampened 
Hertz-Mindlin theory, with or without possible global 
dampening due to the viscosity of the air. We have found 
that there is a dynamic Janssen effect in the sense that, 
at any frequency, approximately one-half of the effective 
mass rests on the bottom of the cavity and one-half on 
the side walls. Thus, these toy models have only quali- 
tative validity. Notwithstanding, we use them to extract 
approximate values for the sound speed and the viscos- 
ity of our granular ensembles. Finally, our simulations 
as well as our experiments on the effects of humidity on 
M{uj) indicate that the dominant microscopic mechanism 
for dampening is at the grain-grain contact level and is 
not due to global viscous dampening. 

We show that this dampening is much larger for high 
humidity systems than for low ones leading us to con- 
clude that the mechanism is related to the viscosity of the 
adsorbed water in the region of the contacts. Such a con- 




FIG. 1: Schematic of measurements: The effective mass, 
M(a;), is determined by simultaneously measuring the accel- 
eration, a{oo), and force, F{oo), on an oscillating cup filled 
with grains. Using the measured M((jj) we show how to pre- 
dict the acoustic properties of a resonance structure having 
the same grain-filled cavity. Here, fg represents a resonance 
frequency within the granular medium, /o represents the reso- 
nance frequency of the unloaded structure, and fn represents 
the resonance frequency of the loaded structure. In general, 
these resonance frequencies are complex- valued. 



elusion accords with the finding in "room-dry" sedimen- 
tary rocks that acoustic attenuation is caused by stress- 
induced diffusion of adsorbed layers of volatile molecules 
pj| and, more recently, by direct measurements of acous- 
tic attenuation in dry and weakly wet granular media 

The organization of the paper is as follows: We de- 
scribe our experimental technique in Section [III There 
are two types of measurements here. In addition to mea- 
suring M{u;) we also measure the resonant frequencies 
and dampening rates of flexural normal modes in a rect- 
angular bar having a grain-fllled cavity in it. We show 
how the measured function M{u;) can be used to com- 
pute the effect of the granular medium on the resonant 
frequency and dampening characteristics of a structure 
in Section lllli This is illustrated schematically in Fig- 
ure ([1]). When the cavity in the bar is fllled with grains 
its own resonance frequency is changed from /o in the 
unloaded state to /r, which is complex- valued, reflect- 
ing the attenuation in the problem. Additionally, the 
effective mass, M{u;)^ generally has its own (complex- 
valued) resonance frequencies, fg. It may happen that 
one or more of these resonances is manifest as subsidiary 
resonances within the bar + grains system, although the 
complex- valued resonance frequency is changed by virtue 
of the bar's compliance. 

Next, in Section [IVl we discuss some rather general 
properties of the dynamic effective mass, considered as a 
causal response function. We specify to a model in which 
we treat the granular medium as a collection of rigid ob- 
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jects interacting via contact forces between them. We 
also specify to simple, continuum mechanical, models. 
Section [Vl is devoted to an analysis of our data in the 
context of these continuum models which are useful for 
understanding features such as the main resonance and 
the high frequency tail. Here, we analyze data on cavities 
filled with simple liquids, as well as data on our granular 
media. In order to get a sense of whether contact damp- 
ening or global dampening is the dominant mechanism in 
granular media we have performed a series of numerical 
simulations, which we report in Section I VII In Section 
I VII I we investigate the effect of differing humidities both 
on M(uj) and on the flex bar resonances. By means of 
these data it becomes clear that the dampening is local, 
due to the viscous adsorbed films of water in the contact 
region. We summarize our conclusions in Section [Villi 



Empty cup and loaded cup data 




II. EXPERIMENTAL PROCEDURE 

A. Effective Mass of Granular Media 

A cylindrical cavity (of diameter 2.54 cm and height 
3.07 cm) excavated in a rigid Al cup is filled with tung- 
sten particles. Each of these particles consists of four or 
five equal-axis particles, of nominal size 100 /im, fused 
together. The large density of tungsten maximizes the 
effects we are studying. The cup is subjected to a ver- 
tical sinusoidal vibration at angular frequency cj; the re- 
sulting acceleration is measured with either one or two 
accelerometers attached to the cup, on the underside to 
one side of center (see below), and the force is measured 
with a force gauge mounted between the shaker and the 
cup. Taking into account the mass of the empty cup, Mc, 
we have 



a{uj) ' 



(1) 



where the effective mass of the granular medium, M{u;) = 
Mi{uj) + iM2{oo)^ is complex- valued, reflecting the par- 
tially in-phase, partially out-of-phase motion of the in- 
dividual grains, relative to the cup motion. As we shall 
see. Ml may be conceptualized as an inertial effect, in 
the conventional sense, except that it depends upon an 
interplay between the masses of the individual grains and 
the stiffnesses of the grain-grain contacts; M2 describes 
the attenuative aspects of the medium. 

Of course, there is no such thing as a perfectly rigid 
cup. We have used a variety of different cups of slightly 
different geometries, with the intention that the appar- 
ent effect mass of the empty cup should be approximately 
frequency-independent. In Figure [2] we show the results 
for one such cup plotted over a fairly wide frequency 
range. We show both the empty cup data (dashed lines) 
and the filled cup data (solid lines). Depending upon 
the way in which the cup is filled with grains, either by 
slightly tapping on the cup (black curve) or by mechani- 
cally compacting the grains with press and plunger (green 



FIG. 2: Measurement results using one specific shaker cup 
over a wide frequency range. The top row represents the mea- 
sured force, the middle row the acceleration, and the bottom 
row the effective mass. The left column represents the mag- 
nitude and the right column the phase relative to the driving 
voltage on the shaker, except that the bottom row shows the 
real, Mi, and the imaginary, M2, parts of M(cj). The dashed 
curves represent the empty cup. The black curves are for a 
cup that was filled by simply tapping it whereas the green 
represents the mechanical compaction protocol, as described 
in the text. The measured effective mass of the empty cup, 
shown as a dashed curve, is essentially a frequency indepen- 
dent constant, which constant is subtracted from the total 
effective mass to yield that of the granular medium. 



curve) we get quite different results, as discussed below. 
We note that the effective mass of the empty cup is es- 
sentially a frequency-independent constant, as expected. 
This result is a consequence of our use of two accelerom- 
eters on either side of the bottom of the cup, and taking 
the average. In either accelerometer there is a visible 
resonance structure around 6 kHz, which seems to be a 
consequence of a "wobble" motion; the cup does not, lit- 
erally, oscillate along a vertical axis. By averaging the 
two accelerometer signals, this wobble motion is effec- 
tively canceled. In this way we have improved on the 
technique reported in Ref. [6]. 

We note from Figure [2] that in the low frequency limit, 
M{(jo) tends to the static mass of the grains. Also, there 
is a relatively sharp resonance peak whose position de- 
pends strongly on the manner in which the grains were 
prepared; there are also subsidiary resonances at higher 
frequencies. Although the data in Figure [2] are for the 
tungsten grains, we have shown previously that these 
general features of M{uo) are present in other granular 
media, such as spherical glass beads or spherical lead 
beads [6]. In this article we focus on the the tungsten 
granules, simply because the effect is maximal for such 
dense particles. 
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B. Resonant Bar 

We consider the resonant frequencies of a rectangular 
bar of stainless steel whose dimensions are L X W X H = 
20.32 cm X 3.81 cm X 3.18 cm. In the center of the top 
surface a cavity is excavated having virtually identical 
dimensions as that in the shaker cup (above). This cavity 
is filled with tungsten granules in the same manner and 
with the same mass as for the shaker. 

We monitor the flexural modes of the system: The bar 
is suspended by wire supports attached at the approxi- 
mate locations of the displacement nodes of the funda- 
mental flex mode of the bar [Table 3.2 of reference [13]]. 
The purpose here is to minimize any additional damp- 
ening in the experiment due to radiation of energy into 
the bar supports. On the top and on the bottom of the 
bar, near each of the ends, we epoxy piezoceramic disks 
of 6.35 mm diameter and 3.18 mm thickness. The two 
top disks are driven in phase with each other while the 
two bottom ones are driven out of phase with the top 
ones. In this manner the four disks act as a bending mo- 
ment on the bar, thus inducing the desired flexural mo- 
tion. On the bottom directly under the cavity we mount 
an accelerometer. The bending moment is driven at a 
constant voltage as the frequency is swept in increments 
which are finely spaced when the Q is large (as for an 
empty cavity) and more coarsely spaced when the Q is 
low, as when the cavity is grain loaded. The output of 
the accelerometer, a{u;)^ shows a characteristic resonance 
peak as we sweep through the resonance. 

We take the power dissipated in the bar to be P{uj) ex 
I a{u;) p /cj^. This quantity is plotted, in arbitrary units, 
in Figure [3l The resonance frequency, /r^ is taken at the 
peak. The quality factor, Q, is taken from the full- width 
at half maximum as indicated in the figure: 

Q = fR/FWHM . (2) 

See pages 23 JJ. of Reference [l3|. In order to check 
that we are indeed seeing a true resonance we examine 
the variation of the phase angle of the acceleration as a 
function of frequency. A true resonance will exhibit a 180 
degree phase change, as is seen in the two examples of 
Figure [3l Moreover, it will also exhibit a 90 degree phase 
change between the two half- max points, which feature 
is also indicated in the figure. 

This technique works well when there is just one peak 
and the background power absorption is small. In the 
cases where the resonance peak occurs in the middle of 
some broad absorption mechanism, or even when there 
are two resonances closely spaced in frequency, we have 
used a different approach. Here, we record the complex- 
valued acceleration data, a{uji), as for the cup, where 
{(jji} is the set of discrete measurement frequencies. The 
desired (complex- valued) resonance frequency is that for 
which the analytic continuation of a{u;) has a pole in 
the complex cj plane. We form the auxiliary function 

g ^= 1/a and we use a rational function technique 



empty bar empty bar 
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FIG. 3: Frequency scan of resonance in bar without (top) and 
with (bottom) tungsten granules. Power / Volt ^ (arb. units) is 
plotted in the left column and the phase of the accelerometer, 
relative to the driving voltage on the benders, is plotted in 
the right. Note that, in this example, there are two distinct 
resonances in the granules loaded bar. The red circles mark 
the frequencies at which the power is equal to half its peak 
value. 

(Bulirsch-Stoer algorithm) to do this analytic continu- 
ation, as described in Reference [14] (See also Appendix 
B). It is relatively simple to search for a zero, uor^ of 
this analytically continued function s.t. g{ujR) = us- 
ing Muller's method [3]. The real- valued resonance fre- 
quency is /r = Real(a;i^)/27r and the Q is given as 

Q = -Real(c^i^)/[2 Imag(c^i^)] > . (3) 

Note that for a well-defined resonance, such as shown in 
Figure [3l our two definitions of Q are equivalent and, 
in fact, we have found that these two techniques yield 
identical results to high accuracy. 

We see from Figure [3] that there is a very significant 
effect due to the addition of the loose grains. When the 
cavity is empty, and the frequency scanned with 0.5 Hz 
increments, we find /o(empty bar) = 3014 Hz ±0.1% and 
Qo(exp) = 1004 ± 15%. When the cavity is filled with 
the tungsten particles, and the bar is scanned in 5 Hz fre- 
quency increments, there is a significant frequency shift 
for the main resonance, /i^(loaded) = 2600 Hz, and in- 
crease in attenuation Qi^ (loaded) =46. It is this lat- 
ter feature that makes granular dampening an attractive 
possibility for reducing unwanted structure-borne sound. 
Moreover, there is a new, additional mode that is not 
present in the unloaded bar. The major thrust of this 
article is that these resonant bar characteristics can be 
computed directly, based on the measured effective mass, 
M(cj); a detailed comparison between theoretical predic- 
tions and experimental results, as a function of humidity. 
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FIG. 4: Photomicrograph of some tungsten granules used in 
the experiment. There is a 125 /xm wire at the right, for 
comparison. 



is presented later in Section IVlI Bl and in Appendix A. 



C. Sample Handling 

A major concern for us is that we need to be able to 
prepare the loose grains in the two cavities in a repro- 
ducible manner. This will allow us to make meaningful 
predictions about the resonant properties of the bar, us- 
ing the measured effective mass in the shaker cup. Gran- 
ular media are notoriously hard to prepare reproducibly 
because they often configure themselves into meta-stable 
states. For spherical grains it is relatively easy to get 
the grains into a state approximating random close pack- 
ing, whose properties are quite reproducible |15|. Our 
granular systems, however, do not lend themselves easily 
to this. This is because each granule consists of several 
sand-grain like particles fused together. A photomicro- 
graph of a few of these granules is shown in Figure [H 
Although the main features of the effective mass are cer- 
tainly reproducible, the details of the structure in M{u;) 
can vary significantly enough from one filling to the next 
that it prevents an accurate theory-experiment compari- 
son in the bar. 

In Figure [5] we show the reproducibility of the effec- 
tive mass for consecutive runs using each of the different 
protocols. In the first row we show the results of filling 
while occasionally tapping to settle the grains. Next, we 
show the result of vibrating the particulate medium at 
a frequency of 1 kHz over a range of acceleration ampli- 
tudes, < a < 30(m/s^) with a 93 gm stainless steel 
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FIG. 5: Effective mass, M{co), measured on subsequent pack- 
ings of the shaker cup using four different loading protocols. 
First Row: simple filling and subjective tapping. Second Row: 
Sinusoidal vibration, aligned cup. Third Row: Sinusoidal vi- 
bration, different cup orientations. Fourth Row: Mechanical 
compaction with commercial press. 



plug resting on the free surface. The plug ensures that 
the grains at the surface experience a static pressure sim- 
ilar to that experienced at the bottom of the particle 
pack. The protocol consists of systematically increasing 
the acceleration amplitude in steps of 5m/s^, holding for 
2 min at each amplitude before increasing to the next 
level. After reaching the maximum acceleration ampli- 
tude {30m/ s'^) the procedure is reversed. Therefore, the 
sample is exposed to 6 different acceleration amplitudes, 
and the protocol consists of 11 total steps. The procedure 
is similar to that employed by Nagel et al. jl5| . 
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The results indicate that vibrating the sample with a 
mass on the free surface yields a fairly consistent M{uj). 
All of the samples exhibit the same features, and these 
features are observed at relatively consistent frequencies. 
The main resonance in the grains is observed to occur 
over a slightly broader frequency band (2.4 - 2.6 kHz) 
than that observed when the sample is tapped. On the 
other hand, compared to that observed after the sample 
is tapped there appears to be slightly better reproducibil- 
ity in the shape and position of the secondary features. 
It is interesting to note that the higher degree of com- 
paction achieved with the vibration protocol shifts the 
main resonance to higher frequency, as compared to that 
observed with the tapping protocol. In Section IV Bl we 
will see that this shift can be interpreted as an increase 
in effective sound speed in the granular medium. 

Unfortunately, while reproducible, the effective mass 
so obtained with this protocol varies with the orientation 
of the cup relative to the shaker on which it is mounted. 
This effect is demonstrated in the third row of Figure 
[5l The effective mass is reproducible for a single ori- 
entation, but it varies with orientation. Evidently the 
wobble effect we discussed earlier is big enough to effect 
the packing of the grains, when we vibrate at these high 
amplitudes. Therefore, the vibration protocol is not op- 
timal for handling the samples prior to a transferability 
experiment because we cannot insure that we can dupli- 
cate the motion of the bar and the cup as we vibrate it 
at high amplitudes during the preparation phase. 

The most successful method for developing a repro- 
ducible loading was to mechanically compact the grains 
after they had been loaded in their respective cavities. 
This protocol consists of using a mechanical testing in- 
strument to impose a sinusoidal stress on the free surface. 
To promote a uniform imposition of the stress over the 
free surface of the grains, we use a stainless steel plunger 
with a rubber pad glued to the bottom surface. First, 
a static stress of 59.2 kPa is imposed on the granular 
medium. Then, a sinusoidal stress is imposed on the sys- 
tem consisting of 200 cycles at a frequency of 0.25 Hz. 
The stress amplitude is systematically varied between 
39.5 kPa and 118.5 kPa, in steps of 39.5 kPa. To prevent 
unloading the system, the static stress is increased by 
39.5 kPa for each equivalent change in stress amplitude. 
After the maximum stress amplitude is achieved, the pro- 
cedure is repeated in reverse. So the first two steps con- 
sist of increasing, and the second two steps decreasing the 
stress amplitude. Overall the system is exposed to 1000 
stress cycles (5 sets of 200) at systematically varied stress 
amplitudes. In order to test the reproducibility of this 
protocol we have repeated the measurement by dumping 
out the grains, repacking the same amount by weight of 
fresh grains using the compaction technique, measuring 
M(cj), dumping out the grains, repacking, etc. For the 
purposes of comparing theory vs. experiment in the bar 
data, we used the mechanical compaction protocol. See 
Section [VHBl 



III. DAMPENING OF STRUCTURE-BORNE 
SOUND AND THE EFFECTIVE MASS OF 
GRANULAR MEDIA 

Subject to the validity of a few simple assumptions, 
it is possible to use the measured effective mass, M{uj), 
to compute its dampening effect on an elastic structure 
having an "identical" grain-filled cavity. Let us suppose 
that there is an acoustic structure which has a reso- 
nant frequency ujq when the cavity is empty. The res- 
onant frequency becomes complex- valued in the presence 
of the granular dampening mechanism: co = coq -\- Au. 
Here, Imag(cj) describes the ring-down rate of the mode 
in the time domain or, equivalent ly, the quality fac- 
tor of the resonance in the frequency domain: Q = 
— Real(cj)/[2 Imag(a;)] > 0. We make the approximation 
that the grains contribute an additional mass-loading 
which is localized at position xi. Thus, the density may 
be written as 

p(x) = po(x) + M{u;)S{^ - xi) (4) 

where po(x) is the point-by-point density of the struc- 
ture when the cavity is empty. This assumption is es- 
sentially a statement that we are considering only those 
normal modes of the structure whose wavelengths are 
much larger than the dimensions of the cavity. We make 
the further assumption that the grain-filled cavity con- 
tributes negligibly to any change in the effective elastic 
moduli of the host structure, i.e. The effective moduli of 
the grains, as seen by the walls of the cavity, are much 
smaller than those of the host material. In Section fV Bl 
we show that this is an excellent approximation for the 
effective masses and structures we are considering here. 

The equations of elasto-dynamics may now be recom- 
puted using this modified, and now frequency dependent 
and complex- valued, density. We demonstrate how to do 
this sort of computation for the specific case of the fun- 
damental flex mode in a rectangular bar, the geometry 
considered in Section [HBl in Appendix A. 

It is most useful, however, to consider the results of the 
much simpler perturbation theory, which is valid to first 
order in M(c<;), considered as a small perturbation. The 
displacement field obeys the usual equations of motion: 

- pXui = [CijkiUk,i]j . (5) 

Here, p and Cijki are the position dependent density and 
elastic constants of the material, Ui{x.) is the i-th com- 
ponent of the displacement field, a comma denotes dif- 
ferentiation w.r.t that coordinate, and summation over 
repeated indices is understood; A = cj^ is the eigen- 
value for the problem. Written in this manner. Equation 
([5]) applies to any position-dependent material constants, 
p(x) and Cijfez(x), including those with step discontinu- 
ities. We consider resonances such that either the dis- 
placement field vanishes on the boundary surface of the 
structure, u \s= 0, or the stress tensor vanishes on that 
surface, CijkiUk^i 15= 0. Moreover, we assume that the 
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elastic constants, Cijkh and the unperturbed density, po, 
are real- valued, thus guaranteeing that cjq, the resonance 
frequency when there are no grains in the cavity, is also 
real- valued. 

Now, if the density is perturbed to take on a new value, 
P = Po + then there is a corresponding change in 
both the eigenvalue A = Aq + A A, and in the eigenvector, 
u = + Au. Substituting into Eq. ([5]) and collecting 
all the first order changes one has 

- AXpou^ - XoPo^Ui - XoApu^ = [CijkiAuk,i]j . (6) 

Multiply Eq. ([6]) by [u^]* (sum on i) and integrate over 
all space. The term J[u^]''[CijkiAuk,i]j dV can be in- 
tegrated by parts; the surface term vanishes because of 
the assumed boundary condition (above). The remain- 
ing volume integral cancels the term — Aq / po[u^YAuidV 
(because uq satisfies the zero-order equation). Thus, 



AA 



Ao/Ap(x)|uQ(x) p dV 
/po(x)|uO(x)P dV 



(7) 



Applied to the case of interest, Eq. (|4]) for which Ap(x) 
M(a;)(5(x — xo), this now reads: 



Acj 



2Mn 



where 



Mo I u°(xi) p 
/^o(x)|uO(x) |2 dV 



(8) 



(9) 



in which u^(x) is the displacement field, and po(x) is 
the density, when the cavity is empty of grains. Mq = 
/ Po(x) dV is the total mass of the structure when there 
are no grains in the cavity. Written in this manner, Ii 
is dimensionless. In order for this sort of theory to be 
valid, the granular media in the two cavities must be 
substantially the same. In Section IVII Bl we make such 
a theory-experiment comparison, but we have had to go 
beyond the simple perturbation theory result for two rea- 
sons: (1) The effective mass, M(a;), is not always small. 
It can, in fact, take on values larger than that of the 
steel bar. Perturbation theory cannot give an accurate 
description in these cases. (2) Some of the modes seen 
in the bar + grains system are basically modes within 
the grains [i.e. poles of M(cj)] modified by the effect of 
the bar. Perturbation theory is silent as to these modes. 
Nonetheless, for the modes which are primarily bar-like. 
Equation (|8|) gives a useful intuitive way to think about 
the effects of the granular medium on the resonance fre- 
quency and Q. Specifically, Mi{uJo) determines the shift 
in the resonance frequency and M2(a;o) determines the 
lowered Q, as is clear from Equation (j3|). 

We note in passing that the effective mass is, in reality, 
a tensor, viz Mij{uj), refiecting the fact that gravity plays 
a major role in establishing stiffness and dampening at 



the contacts. Equation ([5]) has an obvious generaliza- 
tion to the case of a tensorial density and Equation ([8|) 
becomes 



cjQiiQ*(xi)Mi^-(cjo)iiQ(xi) 
" 2/po(x) |uO(x) |2 dV 



(10) 



In this article we are considering only situations in which 
the cavity motion is strictly along the z-axis, parallel to 
the force of gravity. Thus the only relevant component 
of the effective mass being considered here is Mzz{^) 
which we shall henceforth denote without the subscripts, 
but with this understanding that there are other nonzero 
components of the tensor. 



IV. PROPERTIES OF M{uj) 

In this Section we first discuss some general properties 
of M{uj) considered as a causal response function. Next, 
we investigate some properties of the system in which we 
idealize the grains as being rigid particles which interact 
with their neighbors via contact forces, of which there 
may be different kinds. A specific motivation here is to 
analyze the high frequency behavior seen in our own mea- 
surements of M{u;) within this discrete particles context. 
Finally, some of the features we observe in M{u) are sug- 
gestive of a collective motion in which the displacement 
varies slowly from grain to grain. This suggests the pos- 
sible approximate validity of continuum models, two of 
which we present here. 



A. General 

The effective mass, M(c<;), is a causal response func- 
tion in that one may, in principle, apply an arbitrary 
time-based protocol to the acceleration of the cup and 
measure the force induced by this protocol. As such it 
has several general properties, which we summarize here. 
These properties follow on general principles as described 
in e.g. Landau and Lifshitz [l^. M(cj) is the fourier 
transform of a real- valued memory function: 



M(cj) 



X{t)e'^^ dr 



(11) 



Causality considerations specifically restrict the range of 
integration to positive values of r only. That is, in the 
time domain, the force and the acceleration are related 
to each other via the memory function, x(^)- 



Fit) 



f 

Jo 



X{r)a(t - t) dr 



(12) 



Inasmuch as only the past history of the acceleration mat- 
ters in Equation (p!2|) its fourier transform leads to Equa- 
tion (pT|) . If the frequency, uo^ is extended to take on 
complex values, we see that M{uo) is regular everywhere 
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in the upper-half complex plane. We also see from Eq. 
(fTT|) that 



M{-oj*) = M*{u) 



(13) 



where an asterisk signifies complex conjugation. These 
considerations lead immediately to the usual Kramers- 
Kronig relations between the real and the imaginary 
parts: 



Mi(cj) 



M2H = -— P 

TT 



f 

Jo 



xM2{x) 



Mi(x) 



dx 



dx 



(14) 



(15) 



where uo and x take on real values only and P[ ] de- 
notes principal part of the argument. (We are assum- 
ing there are no singularities on the real axis and that 
\\m^^oo M{uo) = 0. See Section HVlB. below.) 

One may just as well consider a time-based protocol 
for an applied force and measure the induced accelera- 
tion of the cup. Therefore 1/M{(jo) is also a causal re- 
sponse function and all the results quoted above apply 
to it as well, except for the Kramers-Kronig relations. 
Thus, M{uo) has no zeroes, poles, branch points, or sin- 
gularities of any kind anywhere in the upper half plane. 
Inasmuch as only the constant functions are analytic and 
bounded everywhere in the complex plane, M{uo) must 
have singularities of some sort in the lower- half plane. 

The power dissipated in the system, averaged over one 
cycle of oscillation, is 



{Real[F]Real[v]) 
{l/2)Real[Fv''] 
(1/2) ljM2{uj) I V p 



(16) 



where v = ia/uo is the velocity of the cup. Inasmuch as 
P is always positive, one has uoM2{ijo) > 0, for real values 
of the frequency. We see that this property is borne out 
by the data in Figures [2] and [5l 



B. Systems of Discrete Particles 

Let us consider a model in which each grain is consid- 
ered to be rigid except for the region near the contacts 
with its neighboring particles. Let be the equilibrium 
position of the center of mass of the i-th particle, whose 
mass is rui , and be its displacement from equilibrium. 
Similarly, Oi is the librational angle of rotation. If two 
neighboring particles translate or rotate such that their 
points of contact would move relative to each other there 
will be a restoring force due to the contact forces. The 
linearized equation of motion for the i-th particle is 



-m^cj^Ui = -Ki^ ■ [ui -^OiX di^ - W] 
Kij • [u_^- + Oj X dji - Ui - Oi X dij] 



(17) 



where dij is the vector from to the point of contact 
with the j-th grain. It is understood that the tensor 
'Kij{= Kj^) is nonzero only for grains actually in con- 
tact with each other. (We assume there is at most one 
contact per pair.) Similarly, d^^^; and K^^ refer to grains 
that are in contact with the surfaces of the cavity, whose 
rigid displacement is W. 

The equation of motion for the angular variables is 

-00% • Oi = -di^ X Ki^ • [ui -^OiX di^ - W] 
+ J2j dij X Kij ■ [uj + dj X dji -Ui-diX dij] 

(18) 

where 1^ is the moment of inertia tensor for the i-th par- 
ticle. 

In the special case that the particles are identical 
spheres we have dij = (1/2) [Xj — X^] and the spring 
constant tensor may be written in terms of normal (N) 
and transverse (T) stiffnesses as 



k^.dijdij ^kjl-dijdij] 



(19) 



where dij is the unit vector and we use dyadic notation. 
Similarly for the contacts with the walls. An example 
here would be Hertz-Mindlin contact forces in which the 
stiffnesses increase with increasing static compression but 
we also consider forces of adhesion, capillarity, etc. 

It is understood that, generally, each of the elements of 
the tensors 'Kij or K^^^ are complex- valued and frequency 
dependent reflecting the microscopic origin of the dissi- 
pation. For example, one may take 



K,j{Lo)=K° 



ico'Bi 



(20) 



in which the second term describes an inter-particle force 
proportional to the difference in velocity of the two 
grains. The tensor By is analogous to a dampening pa- 
rameter. In general Eq. ([2Q|) represents simply the first 
two terms in the Taylor's series expansion of'K{uo). These 
"springs" may have rheological properties of their own. 
For example, if there is an internal degree of freedom 
with relaxation time r it is easy to show that 



K,,{oo)=Kl 



K? 



1 



(21) 



The derivation of Eq. ^ parallels that of Eq. (78.6) 
in Landau and Lifshitz [17]. A Taylor's series expansion 
of Eq. ([ni) for small u has the form of Eq. ([20]) for the 
first two terms. 

Notwithstanding the foregoing remarks there are some 
general conclusions one can draw from Eqs. (pT|) and 
(p!8|) . First, the total force which the cavity exerts on the 
grains is 



F = -E^^^w■[u^^O,Xd,^-W] 



(22) 



where the second equality follows because the interparti- 
cle forces cancel, by Newton's third law, as is clear from 
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Eq. ([T7j). Secondly, one can formally write the effec- 
tive mass in terms of the normal-mode frequencies of the 
system: 

M{uj) = y (23) 

n 

where ujn are the complex-valued frequencies for which 
Eqs. (pT|) and (p!8|) have nontrivial solutions when W 
is set equal to zero. Each matrix represents the 
strength of each resonance. Also, from Eqs. (fT7|) and 
(p!8|) it is clear that when the frequency tends to zero 
one has lim^^o = W and lim^^o = 0. There- 
fore, in this limit one has, from the second of Eq. (|22|) . 
F = — cj^MqW where Mq = total static 

mass of the grains. Therefore, from the definition of the 
effective mass, one has: 

lim M(cj) = Mol , (24) 

which seems obvious. 

The high frequency limit of the effective mass also 
has a simple form. From Eqs. ([T7j) and (p!8|) one has 
limc^^oo = and lim^^oo = 0, because of the dom- 
inance of inertial effects, i.e. The particles do not move 
at all, in this limit. Therefore, in this limit the first of 
Eq. ([22|) implies F = lim^^oo ^iwi!^) ' W which, in 
turn, implies 

The high frequency behavior of M(c(;) is controlled by 
the behavior of those "springs" connecting the particles 
directly to the cavity motion. If, for example, these are 
strictly dampened springs, such as implied by Eq. ([2Q|) . 
then 



lim ^A{uJ) 




which is predominantly imaginary valued in this limit. 
On the other hand, if the individual springs have their 
own rheological behavior, such as that implied by Eq. 
(f2T]) . then 

TZ'OO 

lim M(w) = - . (27) 

Should this latter case hold then comparison with Eq. 
(fT4j) immediately gives an f-sum rule on the absorptive 
part of the effective mass: 

Lolmiig[M{Lo)]dLO=-Y^KZ ■ (28) 

So if those grains that are in contact with the walls of the 
cavity have spring constants which are real- valued in the 
high frequency limit then the integrated absorptive part 
of the effective mass is determined solely by those high 



frequency spring constants, independent of the details of 
the grain-grain contact forces. 

The potential usefulness of the results of this subsec- 
tion is that the high frequency limit of the measured 
M{u;) may be compared against Eqs. ([25|) and its vari- 
ants. In our own experiments we are not able to probe 
frequencies high enough to observe this behavior, as dis- 
cussed in Section IV B( below. 



C. Continuum Models 

On a semi-quantitative level, we may understand the 
general features of our experimental results, such as those 
in Figures [2] and [5l in terms of two over-simplified con- 
tinuum models whose main purpose will be to allow us 
to deduce approximate parameter values from our data. 

Model I: The granular medium is considered to be 
a lossy fiuid, with negligible viscous effects a t the wal ls; 
that is, the viscous skin depth [l^ S = y^2r]/{pcj) is 
small compared to the radius of the cup, a. Here, 77 is 
the viscosity of the fiuid; p is the density. The effective 
mass is simply 

M{u;) =7Ta^ptdin{qL)/q , (29) 

where L is the height of the fiuid column, K = Ko[l—icjr] 
is the (lossy) bulk modulus of the fiuid, and q = uj^J p/K 
is the wave vector in the fiuid. This model presupposes 
that 100% of M{uo) is supported by the bottom surface of 
the cup. For "small enough" values of the damping pa- 
rameter, r, resonance peaks, as seen in Figure [6^, occur 
when qL equals odd multiples of 7r/2 i.e. L equals odd 
multiples of 1/4 wavelength: L = A/4, 3A/4, 5A/4. • • • 
These sharp resonances are examples of the reso- 
nances within the granular medium/liquid. The values 
of Kq and r are chosen to mimic the observed frequency 
position and resonance width, respectively, in the exper- 
iments. There is a second resonance in Figure [6^ around 
4500 Hz, but the width of that resonance is nines times 
as large as the first one, so it is scarcely visible in the 
plot. 

We emphasize that this last feature, the resonance fre- 
quencies being in the ration 1:3:5:7 is an artifact of 
three assumptions in Model I: (a) There is no shear rigid- 
ity or shear viscosity in the material, (b) The sound 
speed is constant throughout the sample, (c) The atten- 
uation parameter, r, is negligibly small. We shall see that 
all of these assumptions are violated under the conditions 
of our experiments on real granular media. 

Model II: The granular medium is considered to be 
a very viscous fiuid, which is infinitely compressible. 
This situation may be approximately correct if the filling 
depth of the cavity is much greater than the diameter so 
that most of M{(jo) is borne by the walls of the cup, little 
by the bottom surface. We find 

M{uo) = 2p7ia^LJi{na)/[najQ{na)\. (30) 
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FIG. 6: (a) Model I, a fluid with a lossy bulk modulus, (b) 
Model II, a highly compressible, highly viscous fluid. 



Here n = ^JipuoJr}^ and Jk{z) is a Bessel function of order 
k. This model gives a broad peak with a slow decay at 
higher frequencies as seen in Figure [6)3, much like that 
seen qualitatively in the granular data above the funda- 
mental resonance. We note, for future use, that the high 
frequency limit of Eq. ([3Q|) is 



(31) 



M{uo) 2iiaL^ irjp/uo . 



This limit holds when the viscous skin depth [16], S = 
^J2r]l puj , is much smaller than the radius of the cup. 
We suppose, then, that the result M{uj) oc y/ijuj holds 
generally, for other geometries, but with different prefac- 
tors. 

Is the high-frequency limit for this continuum model 
consistent with what one might expect for a discrete sys- 
tem of particles, with dampening due to their relative 
motion, namely Eq. ([26|) ? Though one may expect 
viscous-like dampening analogous to that described by 
Eq. (|3T]) . at high enough frequencies the viscous skin 
depth becomes small compared to the inter-particle sep- 
aration and there is a crossover from M oc (1 ^ i)uj~^l^ 
to M (X icj"^. This can be demonstrated explicitly if one 
considers a one-dimensional string of point masses, sepa- 
rated a distance h from each other, each of which experi- 
ences a drag force proportional to the difference between 
its velocity and its neighbor's, viz: F = ^\vj — Vj\. This 
leads to the equation of motion for the ensemble: 

-iijj^{u2^\ — 2uj^U2-\) j = —N^...N (32) 

subject to the boundary condition 

^±(iv+i) = W . (33) 



Equation ([32|) is a simple example of Equation ([2C 
which Bij = 7 if z and j are nearest neighbors and Kij = 
0. It is simple enough to solve for the effective mass 
implied by this toy model. Let 



y± 



2^ 



The dynamic effective mass of one such row is 
2z7 



y 



y 



y 



N 



(34) 



(35) 



If we imagine that there is a sequence of these chains, in 
parallel with each other, all connected to the walls then 
Eq. (|32|) is a discretized version of the linear Navier- 
Stokes equation for which the viscosity is 7^ oc 76. The 
effective mass per unit area of the side wall implied by 
the continuum Navier Stokes equation in this geometry 
is 



Mns{oj)/A = — tan(A^T/2) 



(36) 



where T is the separation between the walls. Equation 
(|36|) has a high-frequency limit analogous to Equation 



(37) 



Taking into account the appropriate normalization of the 
relevant constants one may directly compare Eq. ([35|) 
against Eq. ([36|) . This is done in Figure [71 where we plot 
the imaginary part of the effective mass implied by each 
of the models. They agree with each other over much 
of the frequency range. The peak around 300 Hz occurs 
when the viscous skin depth is approximately equal to 
the wall separation, T. Above that frequency each model 
has a frequency dependence M2 ^ ^ as expected 

from Eq. (|37|) . When, however, the macroscopic viscous 
skin depth, b = ^J2r]/puo^ is approximately equal to the 
inter-particle separation, 6, the discrete model, Eq. ([35]) . 
crosses over to a behavior implied by Eq. ([26]) : M2 oc 
uo~^ . In Figure [7] this crossover is visible around 3 MHz. 

To summarize this Section we may say that there are 
several distinct possible origins to the frequency depen- 
dence of M . Within the context of the discrete model, 
such as embodied in Eqs. ([T7|) and ([T8|) , one expects 
structure near those normal mode frequencies which are 
relatively near the real axis, as implied by Eq. ([23]). 
To the extent that a continuum approximation may be 
relevant for some of the lower lying modes, one may ex- 
pect structure when the (wavelength of sound/ viscous 
skin depth) is comparable to the cavity dimensions, ex- 
amples of which appear in Figures [6^, [6)3, or [71 There 
may be further structure at the frequencies for which the 
continuum theory breaks down, an example of which is 
seen around 3 MHz in Figure [3 And finally, there may 
be structure due to the fact that the springs themselves 
exhibit a non-trivial frequency dependence, an example 
of which is given by Eq. ([2T]) . 
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FIG. 7: Comparison of a continuum versus a discrete theory 
for a one-dimensional model of local contact dampening. The 
discrete model closely follows the continuum result except for 
frequencies high enough such that the macroscopic viscous 
skin depth becomes smaller than the interparticle separation. 



V. EFFECTIVE MACROSCOPIC PROPERTIES 
A. Liquids 

A wide variety of real liquids satisfies the assump- 
tions of Model I. We have measured M{u;) for four sim- 
ple liquids. In Figure [8] we show the results for a com- 
mon, commercially available, fluorocarbon, whose chem- 
ical formula is ^'(CFs (0^2)4)3. By fitting Eq. ([29|) to 
this data, we are able to extract the density, p, and the 
speed of sound: V = ^^/Kojp. These values, measured 
with our effective- mass technique for all four liquids, are 
cross-plotted against those determined by more conven- 
tional means [3,[l^ in Figure [9l As there is a good agre- 
ment, both for density and for sound speed, we conclude 
that our technique for measuring the dynamic effective 
mass M{uj) is an accurate one. 

B. Granular Media 

Encouraged by these results, we naively interpret the 
main resonance in Figures [2] and [5] as being a 1/4 wave- 
length resonance of the compressional sound speed i.e. 
analogous to the resonances predicted by Model I. We 
investigated how this resonance, shifts to higher fre- 
quencies as the filling depth, L, is reduced. Throughout 
the volume of grains in the cup, the sound speed must 
be depth dependent; the gravity- attributed stiffness is 
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FIG. 8: Real part of the effective mass of a fluorocarbon fluid 
compared against the theoretical prediction of Eq. (|29)) . The 
1/4 wavelength resonance is visible around 6 kHz. 
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FIG. 9: Densities and sound speeds of various liquids mea- 
sured with our effective mass technique and crossplot against 
those values determined by more conventional means [Tsl . . 

small at the surface and maximum at the bottom. Nev- 
ertheless, we may estimate an effective sound speed in 
the vertical direction based on this peak frequency and 
on the filling depth of the tungsten particles: 

viL) = ALfg{L). (38) 

Figure [10] shows these estimated speeds as a function 
of filling level, L, for a cavity with diameter D = 2.54 
cm and for cavities of differing diameters, all filled to 
the same level, L = 3.05 cm. For this figure we have 
filled the cavities with two different procedures: (1) We 
vibrate the cavity vertically at 1 kHz with different ac- 
celeration amplitudes, as indicated. (2) We fill simply by 
"tapping gently" on the side of the cup. We note that 
the position of the main resonance in the cup can be very 
different depending upon the filling technique. Nonethe- 
less, these results show the trend of greater speed with 
greater depth, as expected. The values we are report- 
ing ~ 100-300 m/s are of the same order of magnitude 
as those rep orted in the literature using other techniques 

HiiHElHlli- 
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In one case shown in Figure [TO] we first filled the cup 
by the vibration measurement and then lightly tapped 
the side of the cup. This caused the grains to pack less 
tightly and move the main resonance to a much lower 
frequency, and thus a much lower apparent sound speed. 
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FIG. 10: Effective speed of sound in tungsten granules as a 
function of filling depth, L, for constant cavity diameter, D = 
2.54 cm , and as a function of cavity diameter, filled to a con- 
stant depth, L = 3.05 cm. The numbers in parentheses refer 
to the maximum amplitude of vibration during the packing 
protocol. 

Also in Figure [TO] we show our results for the effective 
sound speed in cavities of differing diameters, D, filled 
to a common depth. These data provide evidence of a 
dynamic Janssen effect in the sense that the side walls 
support some fraction of the effective mass. If the side 
walls did not support the effective mass one would expect 
the effective sound speed to be independent of cavity di- 
ameter. We return to this point in our discussions of 
numerical simulations, below. Suffice it to say that the 
results of Figure [TO] rule out the strict applicability of an 
interpretation of the main resonance peak based on Eq. 
([29]) . Model I: The sound speed is a function of depth 
and the material has a shear rigidity. Notwithstanding, 
a sound speed in the range 100-300 m/s implies the elas- 
tic moduli are in the range K ^ 1 — 9 x 10^ Pa, which is 
orders of magnitude smaller than the elastic constants of 
steel. The presence of the granular media in the cavity of 
the resonance bar contributes negligibly to the stiffening 
of that bar, as we have assumed. 

Because at least some of the effective mass is borne by 
the side walls of the cavity, we may estimate the effec- 
tive viscosity of the granular medium based on the high- 
frequency tail of the data, such as in Figures [2] and [5] 
According to Model II, the high frequency tail of M(cj) 
should be proportional to o;"^/^, as in Eq. (|3T|) . when 
the viscous skin depth is small: S{uj) <C a. Roughly 
speaking, this is seen in the data plotted in Figure [TTl 
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FIG. 11: M2(cj) vs. frequency indicating the possible ef- 
fect of a macroscopic viscosity. The simulation based on 
contact dampening is suggestive of a macroscopic viscos- 
ity M2 oc /~^^^ whereas that based on global dampening 
M2 oc is not. The simulation results have been shifted 
vertically. 



Taking into account the prefactors, we conclude that 
the granular medium has an effective acoustic viscos- 
ity ^eff = [3 — 10] X 10^ Poise. This value of rjeff im- 
plies a value of the viscous skin depth S{10kHz) ~ 3mm 
which is both small compared to the cup radius and large 
compared to the grain diameters thus satisfying two of 
the underlying assumptions in Eq. ([3T]) . This value is 
just an order of magnitude estimate as the fluctuations 
around the /"^^^ behavior are quite large and Model II 
is, of course, an oversimplification and an overestimation 
of the effects of the sides of the cavities. Nonetheless, 
if our estimate of the macroscopic viscosity is approxi- 
mately correct that would imply that the crossover from 
M2 (X cj"^/^ to M2 (X uj~^ [i. e. Eq. ([26])] occurs in the 
frequency range [2 — 8] x 10^ Hz. This crossover happens 
when the macroscopic viscous skin depth approximately 
equals the inter-particle separation, (5 = 6, as discussed 
in connection with Figure [T] If this crossover happens at 
all, it would seem to be well outside our experimentally 
accessible measurement range. 

In a similar vein, we may conceptualize the effective 
elastic constants as K{u;) = Kq — ioorj from which we de- 
duce that a crossover from elastic-dominated to viscous- 
dominated behavior occurs at a frequency uOe-v ^ J^o/v 
which is approximately 15 kHz for the granular sys- 
tems we have studied. We may very well be seeing this 
crossover at the higher end of our frequency range in Fig- 
ure [IH 

To summarize this subsection we have found some evi- 
dence that there is a an effective macroscopic viscosity in 
the sense of Eq. ([3T]) being approximately correct. This, 
in turn, suggests that we should not see behavior implied 
by Eq. ([26]) unless the frequency is much higher than 
we are investigating; we don't see that behavior, much 
less that suggested by Eq. ([27|) . Even at our highest 
frequencies the granular system is undergoing collective 
oscillation, albeit a complicated one: It is not the case, 
even at our highest frequency, that only those grains in 
contact with the walls are contributing to M{uj). 
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VI. SIMULATIONS 

The toy models are illuminating as far as they go but to 
obtain a deeper understanding of the dampening mech- 
anism on a microscopic level we have performed molec- 
ular dynamic simulations of M{u;). Here, we consider 
the much simpler system of spherical beads confined in a 
rectangular box. 

In our simulations a static packing at a pre-determined 
pressure is first achieved by methods previously described 
[23, and then we incorporate walls, friction, and the force 
of gravity. The simulations consider the typical Hertz- 
Mindlin contact forces for the normal and tangential 
components, respectively, and the presence of Coulomb 
friction between the particle characterized by a friction 
coefficient fi. See Ref. [23] and references therein for 
a discussion of the underlying physics of these contact 
forces. 

Microscopic dampening is provided by two principal 
mechanisms of dissipation: (a) Local dampening, in 
which the force is proportional to the relative velocity 
between two contacting grains [l^, [HI- Examples here 
would include intrinsic attenuation due to asperity de- 
formation [23], and wetting dynamics within the liquid 
bridges between adjacent grains as they move relative 
to each other [2^, [23] • (b) Global dampening, in which 
there is a presumed rotational and translational drag due 
to a viscous fluid, such as air, which is assumed to move 
with the walls of the cavity. See [25, 28, 29]. Within 
the context of Eqs. (pT|) and ([18]) the global dampening 
approximation is that the particles not in direct contact 
with the walls still experience a drag effect such that 
— — zcjBw, where is the same for each grain. 
The conclusions we draw from our numerical simulations 
are relatively insensitive to the assumed values for the 
dampening parameters, either global or local. We plot 
some typical results in Figure [11] where we indeed see that 
the local dampening hypothesis implies M2 ex o;"^/^ and 
the global dampening is more consistent with M2 ex 

Results are displayed in Figure[T2]for a system in which 
there is assumed to be only local contact dampening. 
The two cases shown correspond to normal and trans- 
verse contact forces (top) and normal forces only (bot- 
tom). The fundamental resonance and the broad high- 
frequency tail are clearly evident. We show, separately, 
the contribution to M{lj) from the bottom as well as from 
each of the side walls. The conclusions we draw from the 
numerical modeling are: 

1. A dynamic Janssen effect reveals that Models I and 
II (Section llV Cp are equally important to an understand- 
ing of the dynamics. Figure [T2k shows the results of the 
effective mass at the bottom and on the walls of the cav- 
ity, as a function of the frequency, for a system with fric- 
tion coefficient /i = 0.5. For all frequencies, we find that 
approximately one-half of the mass is held by the bot- 
tom and the other half by the side walls. In this sense, 
one cannot make the distinction between the simplified 
models. Our depiction of an effective sound speed as a 
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FIG. 12: Computed real part of the effective mass of loose 
granular media, normalized to the static value vs. frequency, 
using molecular dynamics simulations: (Top) /i = 0.5, (Bot- 
tom) /J, — 0. Although both simulations clearly show the 
resonance, in the former approximately one-half of the effec- 
tive mass is supported by the side walls and one-half by the 
bottom. 



function of filling level. Figure [TO] is really just a manner 
of speaking. When the friction is switched off (/i = 0, 
Fig. [T2b). almost all the weight is supported by the bot- 
tom walls of the cavity, as expected, since the effective 
shear modulus becomes negligibly small [l^; therefore 
the Model II effect essentially disappears. (A small com- 
ponent of the weight is held by the walls because they 
are made of glued particles in the simulations.) Not sur- 
prisingly, we still see the resonance peaks as predicted by 
Model I. 

We note that our results for the effective sound speed 
measured in cavities of different diameters but filled to 
the same depth, plotted in Figure [TO] provide an exper- 
imental indication of a dynamic Janssen effect. If the 
side walls did not support the effective mass one would 
expect the effective sound speed to be independent of 
cavity diameter. 

2. Simulations allow us to differentiate between possi- 
ble different microscopic origins of dissipation. We find 
that either global or local dampening can capture the 
main features of the experiments: the main resonance 
peak, as in Model I, and a broad background, as in 
Model II. However, the high-frequency asymptotic be- 
havior for large co is very different for the two mech- 
anisms. Global dampening predicts M{uj) ~ iuj~^, as 
per Eq. ([26]) where Bi^ is nonzero for every particle. 
Roughly speaking this can be seen in the result plotted 
in Fig. (TTJ On the other hand, as long as the viscous skin 
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depth 2r]Qf[ / poo is large compared to the grain size but 
smah compared to the cup radius, then contact dampen- 
ing predicts M{u;) ~ (i/o;)^/^, as in Model 11. This trend 
is seen in the numerical simulations based on contact 
dampening plotted in Fig. [Til Although the simulation 
result is noisy, we may conclude that there is an effective 
viscosity, as is seen in the experimental results. Contact 
dampening can be caused by viscoelasticity of the grains 
themselves or it can be induced by liquid bridges at the 
contact points [2^, [l^, [lO, [3l[ . We are inclined to sup- 
pose it is the latter that dominates in our samples and 
this has motivated us to consider the effects of humidity 
on our results, both for M{u;) and for the resonances in 
the bar. We present our results in the next section. For 
reasons which will become apparent, we need to develop 
a quantitative theory of the bar resonances that goes be- 
yond the perturbation theory, which we also describe in 
the next section 



VII. DETAILED THEORY OF FLEX BAR 
RESONANCES 

The perturbation theory described in Section [TTIl is in- 
formative but in some situations M{ij) takes on very 
large values - larger even than the mass of the bar it- 
self, thus invalidating the assumption that M{uj) is a 
small perturbation. Moreover, there are resonances seen 
in the loaded bar which are primarily resonances within 
the granules themselves; perturbation theory is not able 
to make any prediction about these resonances. For situ- 
ations such as these it is necessary to go beyond the per- 
turbation theory and use a more complete theory. For 
the case of a rectangular bar the flex mode resonances, 
including the effect of the granular medium, can be com- 
puted directly. This technique is described in Appendices 
A and B. Basically, the theory treats the bar as a one di- 
mensional object for which flex waves are described by 
the so-called Timoshenko beam theory [32j. The extra 
compliance near the center of the bar, due to the exis- 
tence of the cavity, is modeled by a single parameter, ^, 
whose value is set by the resonance frequency when the 
cavity is empty. There are no other adjustable param- 
eters. As before, though, we treat the effective mass of 
the grains M{uj) as a localized point density. 

At each (complex- valued) frequency there are two left- 
going and two right-going waves. The amplitudes of these 
four components are determined by the requirement that 
four homogeneous boundary conditions must be satisfied. 
The normal-mode condition is that the determinant of 
these coefficients must vanish; we search, numerically, 
for the complex- valued normal- mode frequencies, {00 r}^ 
at which the determinant vanishes. 

We compare the results of this theory, first against res- 
onance bar data taken when the cavity is filled to different 
depths with various liquids, and then to data taken when 
the cavity is filled with granular media, under conditions 
of differing humidity. 
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FIG. 13: Resonance frequency of steel bar as a function of 
added mass in the cavity, for various liquids, filled to various 
depths. The solid curve is the prediction based on the theory 
described in Appendix A. The dashed curve is the pertur- 
bation theory thereof Eq.(|8]). The two data points indicated 
with squares represent the effective mass of those liquids eval- 
uated at the resonance frequency, using Eq. (|29p . 



A. Liquids 

In Section IV Al we showed how the effective mass mea- 
surements on simple liquids gave results in general agree- 
ment with the predictions of Model I, Eq. ([29|) . It is nat- 
ural to inquire whether the theory described in Appendix 
A can accurately predict the frequency shift of a resonant 
bar whose cavity is filled to varying depths with such a 
simple liquid. In Figure [13] we show the measured results 
of the flexural resonance frequency of the bar when the 
cavity is partially filled with the liquids considered in Fig- 
ure [9l The solid curve represents the predictions of the 
theory described in Appendix A, assuming the effective 
mass within the cavity is frequency independent. Over 
the range of added mass values, it is nearly the same as 
the predictions of the perturbation theory, Eq. ([8|), for 
which we have determined /i ^ 2.1 for the fundamental 
mode in our bar, reflecting the fact that the displacement 
at the center is much larger than the average RMS dis- 
placement. For large enough values of added mass the 
resonance frequency predicted by the full theory asymp- 
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totes to a finite frequency, reflecting the fact that even if 
the center of the bar is pinned, the two arms may stih os- 
cillate freely. The difference between the full theory and 
the perturbation theory is quite large for added masses 
on the order of 500 grams, which is off the scale of Figure 
[9l The data mostly lie on the theoretical curve except for 
the data points greater than 20 gms, which correspond 
to a fluorocarbon fluid. For these data points it is sim- 
ple enough to estimate M(cj), evaluated at the measured 
resonance frequency, using Eq. ([29|) . With this correc- 
tion to the added mass all the data for the simple fluids 
now lie nearly on the theoretical curve, which gives us 
confldence in our approach. 



B. Granular Media - Humidity Effect 
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In order to elucidate the physical origin of the dampen- 
ing mechanism we have undertaken a controlled study of 
the effects of humidity on these systems. Both the shaker 
cup and the resonant bar are fllled with the same amount 
of grains, by weight, of tungsten particles. They are 
packed into their respective cavities using the compaction 
protocol described in Section III CI Both the shaker ap- 
paratus and the resonant steel bar are placed in a her- 
metically sealed glove box. The humidity is controlled 
by placing an open pan of salt-saturated water inside, 
as well. We have used different salts in the water as a 
means of controlling the humidity. We use a low power 
fan to provide a continuous flow of air throughout the 
container. The temperature is held at T = 26.5^ C. The 
glove box sits on a vibration isolation table. The moti- 
vation here is to allow for equilibration of humidity on 
a reasonable time scale and to not allow for extraneous 
vibrations to dislodge the grains in the cavities. 

In Figure [TH we show a comparison of the measured 
and the computed resonance frequencies in the bar as a 
function of elapsed time. We also show how the humidity 
is changing, as we swap one salt-saturated solution for 
another. (The "zero" humidity cases are accomplished 
with a desiccant.) There are several different resonance 
frequencies that were measureable. For each of them 
we have determined the complex-valued resonance fre- 
quency, /i^, using the procedure described at the end of 
Section IIIBi The width of each resonance is indicated 
with the error bars. i.e. What is plotted is Real(/i^) and 
Real(/i^) ± Imag(/i^). Using the simultaneously mea- 
sured M{ijj) in the shaker cup, we are able to compute 
the expected (complex-valued) resonance frequencies in 
the bar, using the theory described in Appendices A and 
B. These computed resonance frequencies are depicted 
in the same manner as the measured ones, in Figure [T4l 
The horizontal dashed line, /o, is the resonance frequency 
in the unloaded bar. Its Q is so high (~ 1000) that the 
width (a few Hz) does not show on the scale of the flgure; 
see Figure [3l top. The other dashed curve, plots the 
position and width of the main resonance frequency seen 
in M(cj). 



FIG. 14: Comparison of measured and predicted resonant 
bar characteristics at various humidities. The "error bars" 
indicate the width of the resonances, as described in the text. 
The filled symbols refer to the dominant mode in the system, 
which is predominantly "bar- like". The horizontal dashed 
line, /o, is the resonance frequency in the unloaded bar. Its 
Q is so high that the width does not show on the scale of the 
figure, fg shows the position and width of the main resonance 
seen in the effective mass data, some of which is show in 
Figure [15] 



A few selected plots of M(cc;) are shown in Figure [T5l 
The main resonance in M{uj) is clearly visible in Figure 
[T5]but, in fact, there is another resonance of smaller am- 
plitude, around 4500 Hz for the 43% and 74% humidity 
cases, which is also visible in Figure [T5l 

There are several interesting features to the data of 
these two Figures. First, from either Figure [Ml or Figure 
[15] we see that the position of the main resonance within 
M(uj) starts at a low (~ 2000 Hz) frequency when the 
humidity is initially zero. It then increases as the humid- 
ity increases to 43% then again to 73%. Concomitantly, 
the width of that resonance also increases; the system 
of loose particles is becoming more dampened with in- 
creasing humidity. Second, there is generally quite good 
agreement between theory and experiment as to the po- 
sition and widths of the several different resonances in 
Figure [TH Third, although the strongest resonance in 
the bar around 3kHz dominates, there are also one or 
two other resonances, which we associate with the dy- 
namics within the granular medium. 

We have noticed that, of the two resonances seen in the 
bar in the region 2000-3500 Hz, one of them has a signif- 
icantly larger amplitude than the other, as is clear from 
the example shown in Figure [3l which corresponds to t 
= 41 hours. This "stronger" resonance is labeled with 
fllled symbols, for both the measured and computed res- 
onances, in the legend of Figure[T4l the others are labeled 
with open symbols. Our interpretation is that the two 
"bare" modes, the resonance in the empty bar at /o, and 
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FIG. 15: Evolution of dynamic effective mass as the humidity 
is changed, as per the humidity data in Figure [Ml 



the resonance in M{uj) at have coupled to each other 
to yield two hybrid modes in the bar + grains system, an 
upper branch and a lower branch. The resonance with 
the larger amplitude corresponds to a mode which is pre- 
dominantly bar-like. When < /o, at early times, the 
real part of the effective mass is negative in the vicinity 
of /o. Although the theoretical calculations were done 
with the complete theory of Appendix A, it is useful to 
consider the qualitative predictions of the perturbation 
theory. According to the predictions of Eq. ([8]), one may 
expect the resonant frequency in the bar should increase 
relative to that of the empty bar due to this negative 
mass loading, and this is just what we see in the up- 
per branch. At slightly later times, when the humidity 
is increased to ~ 70%, fg has increased above /o, with 
the result that the real part of M{u;) is now positive, in 
the region around /q. Now, according to the perturba- 
tion theory, the resonance frequency in the bar should 
decrease and this, too, is just what we see. An analo- 
gous behavior was seen in the data of Kang et al [loj 
who considered the behavior of the resonance frequency 
of a clamped plate as more and more loose beads were 
placed atop it. Our view of their experiments is that for 
low values of mass- loading, the resonance frequency, 
in their grains is higher than the resonance frequency of 
their plate. Thus, the real part of their effective mass is 
positive and the resonance frequency is initially lowered, 
relative to the unloaded value. If the depth of loading 
is large enough, fg drops below the plate resonance fre- 
quency and the real part of M{uo) can become negative. 
This is our explanation of why the observed resonance 
frequency may then increase, relative to the unloaded 
value. 

Although the frequency of the upper branch in Figure 
[T4l is relatively slowly varying, it very abruptly changes 
character, from "bar-like" , at the early times, to "grain- 
like" later. The lower branch shows the opposite behav- 



ior. This aspect of the response of the system will be 
analyzed in detail elsewhere [33| . 

When the system is at 90% humidity, it behaves quite 
differently. (There is an initial drop in humidity which is 
an artifact of our procedure for swapping in and out the 
salt-saturated solutions.) The main resonance in M{ijo) 
initially drops to around 2200 Hz but, more significantly, 
it becomes very much broadened relative to that at the 
lower humidities, cf. Figure [151 We do not understand 
the initial drop in the resonance frequency; conceivably 
it is related to the initial drop in humidity. At this hu- 
midity, we were able to locate only one resonance in the 
experimental bar data and we were also able to locate 
only one in the theoretical computations. As time goes 
on, maintaining the humidity in the glove box constant, 
we see from Figures [14] and [15] that the main resonance in 
M(ijj) increases with time. It increases with time for the 
lower humidities, too, though it is most pronounced for 
90% humidity. Previously [o], we analyzed this aspect of 
the data in the context of an assumed distribution of as- 
perity heights within the contact region. With time more 
and more of the energetically favorable asperity regions 
become filled with condensate, assuming they are smaller 
than the Kelvin radius, Eq. (|39|) . below. Under a fairly 
broad set of assumptions, this lead to the prediction that 
the contact stiffnesses would increase approximately log- 
arithmic in time, which is what was observed. The data 
sets analyzed in Reference [6[ were, however, prepared 
using the vibration protocol described earlier HTCl unlike 
the data in Figures [TH [15] Nonetheless, we believe that 
this is what is happening: At any given level of humidity, 
there is a relatively slow process of equilibration leading 
to a relatively slow increase in the resonance frequency 
within M{uo). 

Finally, when the salt solution is replaced with a des- 
iccant, causing the humidity to drop to essentially zero, 
there is a very dramatic change in the properties of the 
granular medium. According to Figure [15] the main reso- 
nance has shifted up to a very high frequency, ~ 4300 Hz. 
It has become extremely sharp, and extremely strong. 
Concurrently, in the bar there are two distinct reso- 
nances, the one primarily within the grains, around 4500 
Hz, and the other primarily in the bar, around 2800 
Hz. We note that both of these are significantly sharper 
(higher Q) resonances than their counterparts in the hu- 
mid stages. We note that theory and experiment are in 
general agreement as to the level of attenuation of the 
various modes. Specifically, the attenuation of the bar 
resonance (~ 3000 Hz) is much reduced in the dried out 
state relative to that in the humid states. 

We don't really know why we see two resonances in the 
dry cases and one, two, or three in the humid cases, de- 
pending, but this seems to be the situation both for the 
theory and the experiment. In the theory, for 90% hu- 
midity, one may speculate that there is an additional res- 
onance or resonances, presumably related to the apparent 
resonance seen in Figure [T5l but that it /they appear for 
very large values of the dampening rate i.e. well off the 
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FIG. 16: Decay rate, a = — Imagcj^, of the various normal 
modes depicted in Figure [Ml Here, "upper" and "lower" refer 
to the two branches with resonance frequency near 3 kHz. 



real axis. However, if they exist we have been unable to 
locate such resonances with our theoretical machinery no 
matter what value we use for seeding the mode search, 
as described in Appendix A. 

In order to analyze the effectiveness of the granular 
medium in attenuating the modes in the bar, at the var- 
ious humidities, we re-plot part of the data of Figure [TH 
in Figure [T6l Here, a = Imag[cc;i^] represents the decay 
rate of the amplitude, A, of each mode, in the time do- 
main, viz: A (X e~^^. We see that there is a monotonic 
effect of humidity on the decay rate of the main resonance 
within the granular medium, as shown in gold. Such is 
not the case for the dampening coefficient of the mode 
which is predominantly bar-like in the combined system. 
Even though the humidity varies from 43% to 70% in the 
first 80 hours, the dampening rate of this mode exhibits 
only a small, non-systematic, variation, which variation 
is captured both in the experiments and in the theory. 
At t = 80 hours, when the humidity is rased to 90%, 
the measured and the predicted dampening rates of this 
mode both increase enormously. When the salt-saturated 
solution is replaced with a desiccant the measured and 
computed dampening rates drop to very small values, 
much smaller than when the system was originally in the 
dry state. (We note that the grains in the shaker cup 
seem to equilibrate with humidity changes more rapidly 
than do the grains in the bar, but eventually they tend 
to the same state. Perhaps this is a consequence of the 
different flow patterns of air in the vicinity of the two 
cavities in the glove box.) 

Our observed humidity effect is very much analogous 
to that seen in the results of D' Amour et al. [3l| who 
deduced very similar behavior from their measurements 
of the changes in resonance behavior of a quartz plate 
whose top surface was covered with a monolayer of beads. 
They concluded that the in situ drying of the contacts 
between the beads and the quartz causes them to stick 



FIG. 17: AFM image of the surface texture of a tungsten 
granule showing asperities on the 10-20 nm scale. Note the 
differing vertical and horizontal scales. 



over a larger microscopic surface area, causing the asper- 
ities on the surfaces in the contact zone to be partially 
crushed. This, in turn, causes the stiffness of the contact 
"springs" to increase and the dampening to decrease. A 
rough analogy would be when one soaks tissue paper with 
water it is easy to make it stick to a vertical surface; the 
mess is held up by the force of surface tension. If it is 
allowed to dry, the force of adhesion to the wall actually 
increases, even though the water is gone. That is because 
as it dries it also deforms, allowing a large area of solid- 
on- solid contact to develop. This, basically, is what we 
see in our own measurements. An AFM image of the sur- 
face of one of the tungsten granules in Figure [ITl shows 
these asperities on the scale of 10-20 nm. Capillary con- 
densation within these asperities occurs when the Kelvin 
radius, r^, is approximately equal to the asperity size 
where 



RT\n{RH) 



1.04 nm 

ln{RH) 



(39) 



where 7 is the liquid vapor surface tension, R is the gas 
constant, T is the absolute temperature, Vl is the mo- 
lar volume, and RH is the relative humidity The 
numerical value in the second equation is that which is 
appropriate to water at 26.5^ C. Thus, at a relative hu- 
midity of 90% the Kelvin radius is approximately 10 nm, 
which is large enough that the asperities seen in Figure 
[171 have become significantly filled. 

We conclude, therefore, that the dominant mechanism 
of acoustic dampening in granular media is due to the 
adsorbed film of water which exists on the particles, an 
example of contact dampening. In particular, the film 
in the region between two contacting grains undergoes a 
shear deformation due to the relative motion of the two 
grains and thus it gives rise to dampening due to the 
viscosity of the water in the film. The data in Figure [TH 
especially the low attenuation in the second dry state, 
clearly rule out global dampening or intrinsic dampening 
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within the Tungsten as being significant mechanisms for 
dampening. 

It is also clear from Figures [HI and [15] that the ad- 
sorbed films of water at the grain-grain contacts is a 
significant mechanism for stiffness of the inter-granular 
spring constants. However, it is also clear from Figure [TOl 
that something like a Hertz-Mindlin contact forces must 
also be operative: Why else would the effective sounds 
speeds be dependent upon the cavity dimensions? Be- 
cause of the odd shape of our tungsten granules it is dif- 
ficult to analyze the relative contributions to the grain- 
grain stiffness due to each of these mechanisms. For the 
simpler geometry of a random close packing of spherical 
glass beads it is possible to do this estimate, precisely 
because of the data of D' Amour et al. [3l| on monolay- 
ers of beads lying on a vibrating quartz substrate. From 
their published data, their Figure 1 and Equation (5), it 
is simple enough to estimate the strength of the spring 
constants {k in their notation, not to be confused with 
our use of k in our Eq. [30]) between the beads and the 
quartz plate. We find, from their data, 

A^(Single Bead : Exp.) = [1 - 7] x 10"^ N/M , (40) 

depending upon the humidity, in their experiments. 
These values are much larger than what one might expect 
for Hertz contact theory in which the weight of the bead 
provides a static compression and subsequent stiffness of 
the contact. Equations (1-8) of Norris and Johnson [35[ 
in which the normal force, TV, is the weight of a glass 
bead, gives 

A>:(Single Bead : Hertz) = 4.6 x 10^ N/M . (41) 

In doing this estimate we assumed the relevant radius of 
the contact is the radius of the bead itself, 100 /im; had 
we assumed a radius equal to that of a typical asperity on 
the surface of the bead, the computed stiffness would be 
orders of magnitude smaller than Eq. (|4T]) . This estimate 
is further confirmation of the conclusions of D' Amour et 
al [3l| , that the stiffness of the contacts is due to surface 
forces. If, however, one considers the contact stiffness 
of glass beads at a depth of 2.54 cm, compressed by the 
weight of the beads above, one may estimate the expected 
stiffness predicted by Hertz contact theory to be 

/^(Hertz at Depth) = 2 x 10"^ N/M . (42) 

Thus, we may say that the contact stiffness for grains 
near the top surface of a granular-filled cavity - virtually 
any granular-filled cavity - is due to humidity mediated 
surface forces but for those a few centimeters deep Hertz- 
Mindlin contact theory is also important. 

We conclude this subsection by pointing out that the 
measurement technique of D' Amour et al [3l| is very 
much analogous to our measurements on the flex bar. 
In both cases one is looking at changes in the resonance 
frequency of some system due to the, relatively small, 
perturbation of the granular medium. In addition our 



measurements of M{lj) represent a more direct measure- 
ment of the underlying physics of the granular medium. 
In this regard the humidity effect on M{u;) is quite strik- 
ingly apparent, over a wide frequency range as shown in 
Figs. [TH and [151 Although the dampening factor for the 
bar under humid conditions is much greater than that in 
the dry the differences in the main bar resonance are not 
nearly so great as in the effective mass itself. 



VIII. CONCLUSIONS 

We have shown how a measurement of the frequency- 
dependent effective mass of a granular aggregate, M{uj)^ 
allows us to predict, accurately, the effects of a grain-filled 
cavity on the acoustic properties of a resonant structure. 
This fact gives us a more direct access to an investigation 
of the underlying physical mechanisms relevant to the 
dampening effect of granular media on structure-borne 
sound. Crudely speaking, we may think of these systems 
as having an effective speed of sound (small) and an effec- 
tive viscosity (large). The dissipation mechanism occurs 
at the grain-grain contact level as our simulations have 
indicated and our humidity controlled experiments have 
made unavoidably clear. As the humidity is increased 
there is a large increase in the attenuation of the fun- 
damental resonance within the grains, Imag(/^) which 
translates to a non-montonic, but calculable, variation in 
the attenuation of the bar resonance, Imag(/i^). When 
the system is taken to a high level of humidity, and then 
dried, there is a dramatic reduction in attenuation and 
a dramatic increase in stiffness of the grain-grain con- 
tacts, relative to that in the initially dried state. We 
understand this effect in terms of increased solid-on-solid 
contact area at the grain contacts. 
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Appendix A: Treatment of Bar Resonances Using 
Timoshenko Theory 

The very simple one dimensional theory of wave prop- 
agation in a flexing bar [l^ is not accurate enough for 
a quantitative treatment of the effect of a granular effec- 
tive mass on the frequency shift and change of quality 
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factor. Basically this is because the ratio of length to 
thickness of our bar is not large enough (the frequency is 
not low enough). Accordingly we start with the more so- 
phisticated, but still one-dimensional, theory developed 
by Timoshenko [32^ . The equation of motion for the ver- 
tical displacement is 



(A-l) 



where E is the Young's modulus, A = wd is the cross- 
sectional area of the bar in terms of its thickness, and 
depth, d. p is the density of the bar and fi is the shear 
modulus, pi is the moment of inertia of the cross-section 
of the bar relative to its mid-point: / = {l/12)iiP'A for 
a rectangular bar. /c is a shape parameter which we take 
to be equal to (8/9) as is appropriate for a bar of rectan- 
gular cross-section. The simple theory for a flexing bar 
corresponds to keeping only the first two terms on the 
LHS of Eq. (lA^ . 

At each frequency, which is in general complex- valued, 
the solution can be written as a linear combination of 
four linearly independent solutions: 



u = [Asm{qix) + B cos{qix) 



(A-2) 



where qi^2{^) are determined by a direct substitution of 



Eq. (|A-2p into Eq. ()A-ip . Because we are interested 
only in the fundamental mode, which is symmetric w.r.t. 
X = 0, we need only consider the region < x < L/2. 
The coefficients [A, B^C^D] are determined by the re- 
quirement that certain boundary conditions be satisfied 
at the ends on the bars, x = ±i^/2, and at the center, 
X = 0. However, the coefficients A and C are not equal 
to zero, as we shall see, due to the fact that the slope, 
1^, is not continuous at x = 0. 

As we discuss in the main text, we assume the dy- 
namic effective mass is located at the point x = and 
we are treating the bar as effectively a one-dimensional 
object. Taking into account the missing mass of steel in 
the cavity the density is 

p{x) = po + -^-^ -S{x) (A-3) 

where po is the density of steel, Mh = poira'^L is the 
missing mass of steel taken from the cavity, and A is 
the cross-sectional area of the bar. One of the boundary 
conditions is that the net force exerted by the bar on the 
cavity must equal the net mass times the acceleration. 



-[F(0+) -F(0- 
-2F(0+) 



(A-4) 



where F{x) is the force which the bar exerts on the el- 
ement x~^ . Within the context of Timoshenko theory, it 



is specified in terms of specific spatial and temporal time 
derivatives of u{x) [32|]. Thus Eq. (|A-4p gives one of the 
four necessary boundary conditions on the coefficients 
[A,B,C,D]. 

The presence of the cavity at x = makes this region 
of the bar very bendable. The radius of curvature in 
this region is very much smaller than it is away from the 
cavity. We model this effect as if the left and right halves 
of the bar are hinged together such that the bending angle 
between the two halves is proportional to the bending 
moment imposed by the bar: 



9u , >. / 



(A-5) 



where r(x), the bending moment, is also given in terms of 
spatial and temporal derivatives of u{x). The parameter 
^ characterizes the "bendableness" of the cavity region. 
^ = 0^ ^ |o+=0 corresponds to no additional extra 
compliance, whereas ^ ^ oo ^ '^(O^) = corresponds to 
a completely flexible hinge. We determine the numerical 
value of ^ by the requirement that the measured and 
computed resonance frequency of the bar with an empty 
cavity, M = 0, must match. It is the only free parameter 
in the theory. 

The remaining two boundary conditions are that both 
the force and the bending moment vanish at the end of 
the bar: 



and 



F{L/2) = 



r(L/2) = 



(A-6) 



(A-7) 



There are, therefore, four homogeneous linear equa- 
tions in the four unknowns X = [A^ B^C, D] i.e. HijXj = 
and so the determinant of the matrix of coefficients 
must vanish for any non-trivial solution. We search for 
such complex- valued roots, ujr, s.t. det[H{ujR)] = us- 
ing Muller's method [ij]. In general there is a countably 
infinite number of such normal mode solutions. One must 
choose a set of three initial frequencies with which the al- 
gorithm initiates the search. Generally, we have chosen 
the three to be of the form [0.95, 1.00, l.OSjcj^, where the 
seed frequency, cj^, may be complex- valued. The normal 
mode frequency to which the algorithm converges, if any, 
depends upon the value of ujs that is chosen. In the case 
of the data shown in Figure [HI the value of the seed fre- 
quency needed to be very close to the resonances around 
4 kHz, but did not need to be very close to the ones found 
near 3 Hz. 

As a check, we have compared the theory vs. ex- 
periment for a similar steel bar having no cavity at all: 
M = Mh and ^ = 0. The dimensions of this bar are L X 
W X D = 20.32 X 3.05 X 3.81 cm and we used param- 
eters p 7830kg/m^, Vc 5800m/s and Vs 3100m/5 
which we have measured ourselves on the steel samples. 
There are two directions in which this bar can flex; we 
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TABLE I: Comparison of measured and predicted resonance 
frequencies for a bar with no cavity, in Hz. 





Measured 


Full Theory 


Simple Theory 


Flex Horiz. 


3542 


3529 


3793 


Flex Vert. 


4272 


4260 


4741 



measured and we computed the fundamental resonance 
frequency for each orientation. The results are shown in 
Table [H We see that the predictions of the full theory 
agree with the measured values, to within a few tenths 
of a percent (better than the accuracy with which the 
input parameters are known), whereas the predictions of 
the simple flex theory differ by as much as 10%. The 
accuracy of the former approach is more than acceptable 
for our purposes whereas that of the latter is not. 

In order to implement the theory when the cavity is 
filled with a granular medium we need a means to com- 
pute M{u;) for complex- valued a;, though the input data 
is limited to real- valued co. We describe how to do this 
in Appendix B. 

Appendix B: Analytic Continuation of 7kr(u?) for 
Complex Values of u?. 



does propagate. In order to test how well this approach 
actually does the analytic continuation, we consider the 
following simple function: 



F{z) 



1 



sm\z ■ 



(5-i))] 



sm\z ■ 



(7-O.li)] 



}■ 

(B-l) 

We have evaluated this function for 20 different real val- 
ues of and these values are plotted in Figure [TSl 

With these 20 values we construct the rational function 
approximation and plot it for values of z on a line in the 
complex plane parallel to the x-axis: z = x — 0.095i. 

Test of rational function approximation 



Actual function 
- Rational approximation 
Training set on the real axis 




4 6 

Real(z) 



With the shaker apparatus we have determined the dy- 
namic effective mass for a series of real- valued frequen- 
cies: M{ujj) j = 1, • • • , iV. In order to use this informa- 
tion for purposes of computing the resonance frequency 
in the flexing bar, as described in Appendix A, we need to 
be able to evaluate M{u;) for arbitrary complex values of 
CO. This is because the normal mode frequency of the bar 
will be complex- valued, reflecting the attenuation in the 
system. We have chosen to do this analytic continuation 
using the rational function technique, a Bulirsch-Stoer 
algorithm, as described in Reference [IJ. In essence, the 
algorithm approximates the function as a ratio of two 
polynomials - a rational function. The algorithm has the 
property that the rational function passes through each 
input datum exactly. Thus, any noise in the input data 



FIG. 18: Test of rational function approximation. The 
symbols represent the value of the function, Eq. ()B-1|) . 
for values of z on the real axis. These 20 values are used 
to construct a rational function approximation for arbitrary 
complex- values z. It is compared against the original function 
for z = X — 0.095i where we see that there is a near perfect 
overlay of the two functions. 



We also plot Eq. ()B-ip on the same graph where we see 
that the rational function basically overlays the actual 
function, even for values of z near one of the poles. In 
fact, the difference between the two never exceeds 1%. 
This test gives us confidence to use the rational function 
to analytically continue M{u;) into the complex cj-plane. 
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